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Abstract 

We derive closed analytical expressions for the complex Berry phase of an open quan- 
tum system in a state which is a superposition of resonant states and evolves irreversibly 
due to the spontaneous decay of the metastable states. The codimension of an acciden- 
tal degeneracy of resonances and the geometry of the energy hypersurfaces close to a 
crossing of resonances differ significantly from those of bound states. We discuss some 
of the consequences of these differences for the geometric phase factors, such as: Instead 
of a diabolical point singularity there is a continuous closed line of singularities formally 
equivalent to a continuous distribution of 'magnetic' charge on a diabolical circle; differ- 
ent classes of topologically inequivalent non-trivial closed paths in parameter space, the 
topological invariant associated to the sum of the geometric phases, dilations of the wave 
function due to the imaginary part of the Berry phase and others. 

PACS: 03.65.Bz; 25.70.Ef 

1 Introduction. 

During the last ten years the geometric phase factors arising in the adiabatic evolution of 
quantum systems[l] have been the subject of many investigations [2]. Most of the early liter- 
ature was concerned with the geometric phase factors of closed systems driven by Hermitian 
Hamiltonians[3]. More recently there has been substantial interest in the complex geometric 
phase acquired by the eigenstates of open quantum systems. This problem arises naturally 
in connection with various experiments which, by their very essence, require the observation 
of the Berry phase in metastable states ( also called resonant or Gamow states). Dattoli et 
"This work was partially supported by CONACYT (Mexico) under contract No. 4964-E9406 



al.[4] studied the Berry phase in the optical supermode propagation in a free electron laser, 
which is a classical system described by a Schrodinger-like equation with a non-Hermitian 
Hamiltonian. The measurement of the geometric phase in atomic systems with two energy 
levels, one of which at least is metastable, was also described in terms of a non-Hermitian 
Hamiltonian by Miniatura et al.[5]. The validity of the adiabatic approximation for dissipa- 
tive, two level systems driven by non-Hermitian Hamiltonians was examined by Nenciu and 
Rasche[6], and by Kvitsinsky and Putterman[7] who also established that, in this case, the 
Berry phase is complex. A higher order adiabatic approximation for two level non-Hermitian 
Hamiltonians was proposed by CP. Sun[8], who also showed that the holonomy structure 
associated to the Berry phase factor of the non-Hermitian case is the non-unitary generaliza- 
tion of the holonomy structure of the Hermitian case. In a previous paper, we gave closed 
analytical expressions for the geometric phase of true resonant states [9] defined as energy 
eigenstates of a Hermitian Hamiltonian which satisfy purely outgoing wave boundary condi- 
tions at infinity [10], and pointed out some of the mathematically interesting and physically 
relevant properties resulting from the non-selfadjointness of the problem. In this connection, 
later, we showed that the codimension of the accidental degeneracy of n resonances and the 
topological structure of the energy surfaces close to a crossing of resonances differ significantly 
from those of bound statesfll]. By means of a numerical analysis of the experimental data 
on the 2+ doublet of resonances with T = 0, 1 in the energy spectrum of the nucleus 8 Be at 
about 16.6, and 17MeV, we showed in a realistic example, that a true crossing of resonances 
mixed by a Hermitian interaction may be brought about by the variation of only two real 
independent parameters [12]. H-H. Lauber, P. Weidenhammer and D. Dubbers[13] measured 
the geometric phase associated to a triple degeneracy of resonances in a flat, non-symmetric, 
triangular microwave resonator. 

In this paper we explore further the properties of the Berry phase of resonant states of a 
Hermitian Hamiltonian with non-selfadjoint boundary conditions. The paper is organized as 
follows: In section 2, we consider the time evolution of a quantum system in a state which 
is a superposition of bound and unbound states evolving in the presence of an external field 
of force which changes slowly with time. By means of an expansion of the wave function 
in terms of bound and resonant states of the unperturbed system, we associate a complex 
non-Hermitian matrix to the time evolution of the bound and resonant states evolving under 
the action of the time dependent perturbation. Accidental degeneracy of two resonances 
and the properties of the energy hypersurfaces at the crossing of resonances are discussed 
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in section 3. In section 4 we derive general expressions for the Berry phase of resonant 
states and discuss some of its properties. In section 5 we briefly discuss the validity of the 
adiabatic approximation. The paper ends in section 6 with a short summary of results and 
some conclusions. 

2 Adiabatic Mixing of Resonant States 

Let us consider the adiabatic time evolution of a quantum system in a state which is a super- 
position of unstable states moving in some strong external field of force which changes slowly 
with time. In order to have some concrete example in mind, although a very hypothetical 
one, we may think of an 8 Be nucleus which has only unstable energy eigenstates moving in 
the field of a heavy doubly magic nucleus, like 208 Pb, in a peripheral collision in which the 
distance between the two nuclei is never much smaller than the sum of the nuclear radii. In 
a semiclassical treatment of the collision, when the centers of the nuclei move along given 
classical trajectories, the parameters in the nucleus-nucleus interaction change with time[15]. 
Another example could be a highly excited Hydrogen atom in strong external, crossed electric 
and magnetic fields. This system has no bound states, it has only resonances [16]; when the 
external fields change slowly with time, the parameters in the atom-field interaction change 
with time. 

The evolution of the unstable quantum system under the influence of the external per- 
turbation is governed by the time dependent Schrodinger equation 

,n- = «. (i) 

The Hamiltonian H is the sum of the time-independent Hamiltonian Ho describing the evolu- 
tion of the unperturbed unstable system plus a perturbation term Hi parametrized in terms 
of some set of external parameters{Ci, C2....CV}- 

The energy eigenfunctions of the unperturbed Hamiltomian are the solutions of the equa- 
tion 

H cp m (£i) = £ m ip m {£i) (2) 

plus the appropriate boundary conditions. 

Assuming that in the absence of perturbation the unstable nucleus decays spontaneously 
in two stable nuclear clusters, the unperturbed energy eigenfunctions may be written as 
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cluster model wave functions [17], 

= A{M^i)MQ uMrAB) yj I (r AB )) , (3) 

I r AB ) 

where <pA and <Pb are the wave functions of the clusters A and B, uab^^ab) is the radial 
part of the wave function of the relative motion of the two clusters yj 1 {tab) is a spherical 
harmonic, and A is the antisymmetrizer. In our example, ip m would be the eigenfuction of 
a state of 8 Be which decays spontaneously in two A He clusters. In this case, UAB,m{fAB) is 
a Gamow function which vanishes at the origen and behaves as a purely outgoing wave for 
large values of the relative distance tab, 

1 / du AB ,m{r AB )\ ... 
7 V ~a — > lkm ' W 

u AB,m{ r AB ) \ dr J r AB ^oo 

k m is the complex wave number. This boundary condition is not self-adjoint, in consequence 
the Gamow states are not orthogonal in the usual sense, and the energy eigenvalues are 
complex, with Re£ m > and Im£ m < 0. Nevertheless, bound and resonant states form a 
biorthonormal set with their adjoints, which may be extended by a continuum of suitably 
chosen scattering states of complex wave number k to form a complete set in which any 
square integrable function may be expanded. The rules of normalization, orthogonality and 
completeness satisfied by the eigenstates of the Schrodinger equation with local and non-local 
potentials belonging to complex eigenvalues with ImE n < 0, (Gamow or resonant states) are 
given in A. Mondragon, E. Hernandez and J.M. Velazquez- Arcos [10]. A clear exposition of 
the properties of resonant states may be found in the textbook by A. B6hm[14]. 

Since we are interested in the time evolution of the system in a state which is a super- 
position of unstable states, we make an expansion of the wave function * in terms of bound 
and resonant states of Ho, 

* = £o m (tKKi) + / b(k;t)^+\k;0dk (5) 

In general, the index m runs over bound and resonant states. The scattering states 
(/?( + ) (k; £) of complex wave number k are defined by analytic continuation [10]. The integration 
contour C in the complex wave number plane is a straight line with slope -1, that goes through 
the origin. In this way we separate out of the continuum the resonance contribution to 
which should be dominant for energies close to the real part of the complex energy eigenvalues. 

Substitution of (5) in (1) gives the set of coupled equations 
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^# = -i£ m a m (t) - | En < VmltfllVn > «n(t) 

-h/c < >b(k;t)dk 
and a similar expression for ^^'^ , 

When the interactions are time reversal invariant, the dual of the complex Gamow function 
UAB,m{rAB) is the same function [10]. Since the bound state cluster wave functions <\>a and 
4>b are real, the dual of the resonant state wave function ip m (£i) is the same complex function 
fmi^i)- When the interactions are not time reversal invariant, the Gamow function and its 
dual are not the same function. We will use the notation |</?m(£i) >for the Gamow function, 
and < (p m (€i)\ for its dual. Hence, the matrix element of the perturbation term Hi taken 
between bound or resonant states of the unperturbed system is given by 

((PmlH^tpn) =J J {Vm^i)\H 1 \ip n (^))d^ 1 ■ ■ ■ -d^i (7) 
The matrix H with matrix elements 



+ (ip m {H^ipn) (8) 

is, in general, non-symmetric, complex and non-Hermitian. When the forces acting on the 
system are time reversal invariant, H is symmetric, but non-Hermitian. The dependence of 
Hi on the external parameters is inherited by H which is also parametrized in terms of the 
same set {Ci, C2, ■■■■C n } of external parameters. 

In the following we will be concerned with the adiabatic evolution of a system with two 
resonant states very close in energy which are strongly mixed by the Hermitian interaction 
Hi. In our hypothetical example, these could be the two 2 + states of 8 Be with T = 0, 1 at 
Ei = 16.622MeV and E2 = YJ.QlMeV . The contribution of the non-resonant background 
integral over the continuum of scattering functions will play no part in the following discus- 
sion. Therefore, to ease the notation, we will disregard the contribution of the background 
scattering functions. In this approximation, 



m 

Then, the set of eqs.(6) reduces to 



WmiO > ■ (9) 



a m (t) = - T^^-m, n a n {t) (10) 



n 
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where a m (t) is the time derivative of a n (t) . 

The complex energy eigenvalues associated to the resonant states in the presence of the 
perturbation will arise as the eigenvalues of the non-Hermitian matrix H associated, through 
eq.(10), with the decaying process. 

3 Accidental Degeneracy of Resonances 

Non-trivial geometric phase factors of the energy eigenvectors or eigenfunctions are related 
to the occurrence of accidental degeneracies of the corresponding eigenvalues [1, 19]. In the 
previous section we associated a complex nxn, non-Hermitian matrix H to an open quantum 
system with n resonances and resonant states. In the absence of symmetry, degeneracies are 
called accidental for lack of an obvious reason to explain why two energy eigenvalues, Ei 
and E2, of H should coincide. However, if the matrix H is embedded in a population of 

complex non-Hermitian matrices {H(Ci, C2, , CW)} smoothly parametrized by N external 

parameters (Ci, C2, CW), degeneracy in the absence of symmetry is a geometric property of 
the hypersurfaces representing the real or complex eigenvalues of H in a (iV + 2)— dimensional 

Euclidean space with cartesian coordinates (Ci, C2, , Cn, ReE, JmE). In contrast with the 

case of Hermitian matrices, complex non-Hermitian matrices with repeated eigenvalues can 
not always be brought to diagonal form by a similarity transformation. This feature leads 
to a richer variety of possibilities. Thus, it will be convenient to examine first the accidental 
degeneracy of two resonances and the topology of the energy surfaces close to a crossing of 
resonances in parameter space. 

3.1 Degeneracy of two resonances. 

Since we are interested in the possibility of an accidental degeneracy in a system with two 
resonant states strongly mixed by a Hermitian interaction, all other bound or resonant eigen- 
states being non-degenerate, we may suppose that we already know the correct eigenvectors 
of H for all the real and complex eigenenergies E s , except for the two the crossing of which 
we want to investigate. Using for this two states two vectors which are not eigenvectors but 
which are orthogonal to each other and to all other eigenvectors, we obtain a complete basis 
to represent H. In this basis H will be diagonal except for the elements H12 and H21. The 
diagonal elements Hn and H22 will, in general, be non-vanishing and different from each 
other. There is no loss of generality in this supposition, since any complex matrix H may be 
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brought to a Jordan canonical form by means of a similarity transformation[18]. When the 
eigenvalues are equal, H 2X 2 is either diagonal or equivalent to a Jordan block of rank two. 
Hence, in the following, we need consider only the conditions for degeneracy of the submatrix 

H2x2- 

It will be convenient to measure the resonance energies, Eiand E2, from the centroid, £, 
of the diagonal terms in H 2x2 , then, 

H 2x2 = £l + H, (11) 

the traceless 2x2 matrix H may be written in terms of the Pauli-matrix valued vector 
a = (ai,a 2 , 03) as 

H = (12) 



2 

where R and F are real vectors with cartesian components (X\, X 2 , X^&nd (Ti, F 2 , IT^given 
by 

X t -i l -T l = l -[Tr.{7ia, i )}, (13) 

when the forces acting on the system are time reversal invariant, X 2 and r 2 vanish. 

In the absence of more specific information about the external parameters Cj, we will 
parametrize H in terms of R and T , according to (11) and (12). It will be asumed that 
Xi = Xi(C±, C 2 , ....Cn) are bijective and define a homeomorphism of the manifold of the 
external parameters onto the manifold of the complex, non-Hermitian 2x2 matrices H 2x2 . 

From (12), the eigenvalues of TL are given by 



e = T\l[R-i 1 J) 2 (14) 
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and the corresponding eigenvalues of H are 



Ei >2 =£=Fe- (15) 

Then, Ei and E 2 coincide when e vanishes. Hence, the condition for accidental degeneracy 
of the two interfering resonances may be written as 



R d -i-T d ) =0 (16) 
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Since real and imaginary parts should vanish, we get the pair of equations 



R 2 d -\r 2 d = o, 



(17) 



and 



R d ■ r d = o. 



(18) 



These equations admit two kinds of solutions corresponding to H2 X 2 being or not being 
diagonal at the degeneracy: 

i) When both R d and T d vanish, eqs. (17) and (18) define a point in parameter space, H 
is diagonal at the degeneracy and the submatrix H2 X 2 has two cycles of eigenvectors of lenght 
one. The two complex eigenvalues of H which become degenerate migrate to the real axis 
where they fuse into one real positive energy eigenvalue embedded in the continuum [21]. Since 
in this case all the cartesian components of R d and T d should vanish, the minimum number 
of external parameters that should be varied to produce a degeneracy of two resonances to 
form a bound state embedded in the continuum is four or six depending on the quantum 
system being or not being time reversal invariant. 

ii) In the second case, when the degeneracy conditions (17) and (18) are satisfied for non- 
vanishing Rd and T^, these equations define a circle in parameter space. In this case TCd does 
not vanish at the degeneracy and the corresponding is non-diagonal at the degeneray, 
its first Jordan block H2x2 is of rank two and has one cycle of generalized eigenvectors of 
lenght two, all other Jordan blocks have one cycle of lenght one. In this case the two complex 
eigenvalues of H which become degenerate fuse into one repeated complex eigenvalue £. Since, 
the two linearly independent conditions (17) and (18) should be satisfied for non-vanishing 
values of R d and T d , at least two real, linearly independent parameters should be varied 
to produce a rank two degeneracy of resonances. Hence, the codimension of a second rank 
degeneracy of resonances is two, independently of the time reversal invariance character of 
the interaction. 

If we had considered the possibility of H having more than two equal eigenvalues, we 
could have had degeneracies of higher rankfll, 18]. Triple degeneracies of resonances in the 
absence of symmetry have already been observed by Dubbers[13] in a metallic cavity excited 
by a wave guide. Explicit expressions for the codimension of a resonance degeneracy of v 
resonance eigenenergies for v > 2, have been given by Hernandez and Mondragonfll]. 
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3.2 Energy surfaces in parameter space. 

Close to a degeneracy, the energy difference, e, is given by eq. (14). From this expression, 
the real and imaginary parts of e are 



and 




These equations define two hypersurfaces in parameter space. We are interested in the 
shape of the energy surfaces in the neighbourhood of a crossing resulting from the accidental 
degeneracy of two resonant states. 

First, we shall consider the accidental degeneracy of two resonances leading to a rank 
one degenerate Hamiltonian matrix. In this case, the conditions for accidental degeneracy, 
eqs.(17) and (18), are only satisfied for vanishing R and V. As shown in section 3.1, the 
submatrix H2 X 2 should have at least four free real linearly independent parameters to bring 
about a degeneracy of this type. There are only three independent parameters in R. Hence, 
f cannot be a fixed vector. Since R and T should vary independently of each other, it will be 
convenient to understand eqs.(19) and (20) as defining the energy hypersurfaces in an eight 
dimensional Euclidean space, £%, with Cartesian coordinates {X, Y, Z, u, v, w, Ree, Ime} . The 
coordinates (X, Y, Z) are the Cartesian components of R, while (u, v, w) are those of T. In 
this representation the crossing takes place at the origin of coordinates, that is, where R = 
and f = 0. 

When the Hermitian and anti-Hermitian parts of TC commute the problem simplifies. 
From 



R-a,f-a =t(flxf)-a, (21) 

it follows that, when R ■ T and f • a commute, RxF vanishes and R ■ F is equal to i?r. Then, 
the equations of the energy, (19) and (20), take the simple form 



Ree± =±\R\ = ±^X 2 + Y 2 + Z 2 



(22) 



and 



l me± = T - |r| = T-Vu 2 + v 2 + w 2 (23) 

In this particularly simple case the two hypersurfaces which represent the real and imag- 
inary parts of the energy are double cones lying in orthogonal subspaces with the apices at 
the origin. 

In general R • a and F ■ a do not commute and R x T may be written as RT sin 9 with 
sin6> / 0. To study the behaviour of the energy hypersurfaces close to the crossing point, we 
will approach the origin of coordinates keeping R ■ T/RT = cos 6 fixed and let R and T go to 
zero in such a way that the ratio R/T is constant. Expanding the right hand sides of (19) 
and (20) in powers of {^y^J when R < ^T, or in powers of (fjj) when R > ^T, and keeping 
only the terms of lowest order, we obtain 

Ree± ~ ±a \R\ = ±aVX 2 + Y 2 + Z 2 (24) 

Ime± ~ T /3- \T\ = tP^Vu 2 + v 2 + w 2 (25) 

where the factors a, (3 are R ■ T and 1, respectively, if R < ^T; when R = ^T, a and (3 are 
both equal to \J \R.f ; finally, when R > ^T, a and (3 are 1 and R ■ f respectively. 

Therefore, close to a crossing of rank one, the hypersurfaces representing the real and 
imaginary parts of the energy are two double cones lying in orthogonal subspaces with their 
apices at the same point, which, for this reason may be called a double diabolical point. 

Now, let us consider the shape of the energy hypersurfaces close to a degeneracy of two 
resonances of rank two. In this case, the conditions for accidental degeneracy, eqs.(17) and 
(18), are satisfied for non-vanishing values of R and T and, as shown in 3.1, we need at 
least two free real parameters in H to bring about the degeneracy. For definiteness, we will 
keep the anti-Hermitian part of the Hamiltonian matrix fixed and let the parameters of the 
Hermitian part of H 2x2 vary. Then, T ■ a is a constant matrix, T is a fixed vector and R may 
vary. To simplify the notation, it is convenient to choose the OZ axis aligned with T. This 
may be acomplished by means of a similarity transformation of TL which diagonalizes T • a. 

In the case under consideration, equations (19) and (20) define two hypersurfaces in a 
five dimensional Euclidean space, £5, with Cartesian coordinates {X, Y, Z, Ree, line}. The 
hypersurfaces representing Ree and Ime are in orthogonal subspaces and touch each other 
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only when both Ree and Ime vanish. Therefore, the set of points in the energy hypersur- 
faces corresponding to a degeneracy are all in the subspace £3 with Cartesian coordinates 
{X,Y,Z, 0,0}. 

The second degeneracy condition, eq.(18), requires Rd to be orthogonal to T, and defines 
a plane II in S3 and a hyperplane II in £5 with cartesian coordinates {X, Y, 0, Ree, Ime}. The 
first degeneracy condition, eq.(17), is the equation of a circle of radius in the IT-plane, 
which we will call the diabolical circle. 

In order to have an idea of the shape of the energy hypersurfaces in the neighbourhood of 
the crossing, let us consider the surface resulting from the intersection of the e— hypersurface 
and the ft— hyperplane. The equations of this surface are obtained putting R ■ T = in (19) 
and (20), then 



Ree± 



R 2 - 



+ R 2 - 




(26) 



and 



Ime± = =F 



Hence, when R 2 > jT 2 , 



r 2 - ir 2 




(27) 



(28) 



and 



Ime± = 0, (29) 

the energy difference is purely real. Eq.(28) defines a hyperbolic cone of circular cross section 
in parameter space, 



X 2 + Y 2 - (Ree) 2 = -T 2 
4 



(30) 



As shown in Fig. 1, the diabolical circle is at the narrowest cross section or "waist" of 
the cone where Ree vanishes and the two complex energy eigenvalues are equal. 
Similarly, when R 2 < \T 2 , 



Ree± = 



(31) 
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and 



Ime± = =f ( 



1, 



R 2 J (32) 
now, the energy difference e is purely imaginary. This is the equation of a sphere of radius \ 

r, 



X 2 + Y 2 + (/me) 



-n2 



-r 2 



(33) 



The diabolical circle is at the equator of the sphere, where the energy difference e vanishes, 
and the two hypersurfaces touch each other, see Fig 2. 

Off degeneracy [r? / |T 2 ) , the matrix H that mixes the two interfering resonances is 



n 



(34) 



e -v 

where £ and rj are short hand for X — iY and Z — i^T respectively. H has two right and two 
left eigenvectors and may be diagonalized by a similarity transformation 



K~ X HK 



-e 
e 



\ 



(35) 



where 



K 



(36) 



When the conditions for a degeneracy (i? 2 = |r 2 ) is satisfied, det if vanishes and the 
matrix K^ 1 no longer exists. Therefore, at the degeneracy Hd cannot be diagonalized by 
means of a similarity transformation. In this case, TL takes the form 



2 



1 -i e-** ^ 



V 



,i<t> 



(37) 



/ 



This matrix has only one right eigenvector \£d) and one left eigenvector {£d\, belonging to 
the eigenvalue = 0. It also has one generalized right eigenvector \£d) and one generalized 
left eigenvector (£d\ belonging to the same eigenvalue = 0, and such that 



H d \£d) = 0\£ d ) + \£ d ) 
12 



(38) 



and a similar expression for (£ d \. In consequence, the eigenvectors of 7i d and do not form 
a complete basis. We may add the generalized eigenvector to the set of eigenvectors to have 
a complete basis. Then, any vector may be expanded in this basis. However, one should keep 
in mind that, although \£ d ) and are orthogonal to all left eigenvectors belonging to all the 
other non-degenerate eigenvalues, the orthonormality rules for the degenerate eigenvectors 
are 

(£ d \£ d ) = (£ d \£ d ) = (39) 

and 

(£ d \£ d ) = (£ d \£ d ) = 1 (40) 
At degeneracy, that is, on the diabolical circle, 7i d is equivalent to a Jordan block of rank 

two 



o -r 

M~ 1 H d M= I 2 i (41) 







where 



( e-** \ 
M= 4 , (42) 

xe 2 e 2 J 

Therefore, although H and H are continuous functions of the parameters (X, Y, Z) for 
all values of X , Y and Z, both the Jordan normal forms and the similarity transformations 
leading to them are discontinuous functions of X , Y and Z on all points on the diabolical 
circle. 

In brief, we have shown that, in the case of a resonance degeneracy of rank two, close 
to the crossing, the energy surface has two pieces which lie in orthogonal subspaces. The 
surface which represents the real part of the energy has the shape of an open sandglass 
or diabolo, with its waist at the diabolical circle. The surface representing the imaginary 
part of the energy is a sphere. The two surfaces are embedded in orthogonal subspaces but 
touch each other at all points on the diabolical circle. To the points on the diabolical circle 
correspond degenerate matrices with one Jordan block of rank two, and to the points off the 
diabolical circle correspond matrices with simple, that is, non-degenerate eigenvalues. It is 
convenient to recall that when there are multiple eigenvalues, the reduction of a matrix to the 
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Jordan normal form is not a stable operation. Indeed, in the presence of multiple eigenvalues 
an arbitrarily small change in the matrix may change the Jordan normal form completely. 
However, when we are dealing with a family of matrices depending on parameters, multiple 
eigenvalues are unremovable by a small perturbation. In this latter case, we can reduce every 
individual matrix of the family to a Jordan normal form, but, both this normal form and 
the transformation leading to it depend discontinuously on the parameters. Therefore, the 
diabolical circle is a continuous line of singularities of the family of matrices, unremovable by 
a small perturbation. 

4 Berry Phase of a Resonant State 
4.1 Geometric Phase of a Resonant State. 

After having examined the topology of the energy surfaces close to a crossing of resonances, 
let us go back to the expansion of the wave function ^ of the perturbed system in terms of 
bound and resonant energy eigenstates of the unperturbed Hamiltonian, 

* = El^(&)>a m (t)+ / \<P W {k;ti))b(k,t)dk, (43) 

the summation runs over all bound and resonant states of i?o. As before, in the following 
we will disregard the non-resonating background due to the integral over the continuum of 
scattering wave functions of complex wave number. In this approximation 

H9 =E{l^mte))E H --«(^W) (44) 

m n 

When the collective parameters change slowly with time, the perturbation term H\ (d(t)) 
is an implicit function of time. In this case the complex non-Hermitian matrix H is also an 
implicit function of time 

H m ,„ = eW>8 m ,n + <¥>m|#i(a(i))K> (45) 

For values of the external parameters off the diabolical circle H has no repeated eigen- 
values, in consequence, it may be brought to diagonal form by means of a similarity trans- 
formation 

K _1 HK = E, (46) 
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where E is the diagonal matrix of the energy eigenvalues. The columns in the matrix K are 
the instantaneous right eigenvectors of H given by 

H(C i (t))\j> {s) (t))=fi s (t)\4> ( i ) (t)) (47) 

In an obvious notation 

K=(|^)),|0( 2 )), (48) 

The rows in K _1 are the corresponding left eigenvectors of H, properly normalized. 

(<f> (i) \<t> {j) ) = 5 ij . (49) 

With the help of K we obtain the adiabatic basis, {\ip s (£; Ci(t)))} , of instantaneous bound 
and resonant states of the complete Hamiltonian H, 

Ci (*))> = ]T Wm^i))^m,s{Ci{t)) (50) 
m 

and their adjoints 

{M&Ci(t)\ = Y( K ~\t))sn(<Pn&)\ (51) 
n 

We may now write the expansion of ^ as an expansion in instantaneous energy eigen- 
functions {!&(&;<?,(*)))} of H 

^ = Y\^^i,Ci(t)))a s (t) (52) 

s 

where 

a s (t) = J2(K-\t)) sn a n (t) (53) 

n 

Similarly, the expansion of becomes 

h* = Y \<Ps (6; Q(t)))4(t)a s (t) (54) 

Substitution of (52) and (54) in the time dependent Schrodinger equation(l) gives the set 
of coupled equations 
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>kl J ' ' + m^R&n) ■ ^a m (t) = -i£ s (t)a s (t) (55) 



dt ^/ r "' ™ dt 

m=l 



It will be assumed that the non-adiabatic transition amplitudes are very small 

<< 1, m ^ s, (56) 



1 

la* 



dR 

{<Ps\VR<Pm) ■ -^0"m 



Then, we can make the approximation 

1 da s .a ,.s , . , _ . , dR 
H 

Integrating both sides, we obtain 



--^W-(^|V^).-, (57) 



o a (t) = e " ^o 4( *' )d *'e i7s a s (0) (58) 

where the first factor is the complex dynamical phase, whereas the second one is the complex 
Berry phase given by 

7s =i f{0 s \V R $ s ) -dR, (59) 

J c 

in this expression c is the path traced by the system in parameter space when t' goes from 
to to t. 

The right hand side of (59) may be written as a surface integral with the help of Stokes 
theorem 



= * J2 / / {<fs\^ R<Pm) X (<Pm\VR<Ps) ' (60) 

where £ is a surface bounded by the curve c. This expression may be written in a more 
convenient form by means of the identity 

(<p s \V R <p m ) = a _ a (&|V 'RHi\ip m ) , (61) 

then, 

= . f f {(p s \V R Hi\(p m ) x {<p m \V RHi\<p s ) ■ dE 

7s_l i;^^ (4-4,) 2 _ ' 

provided the surface does not cross the diabolical circle where the denominator vanishes. 
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4.2 Computation of the geometric phase. 

Explicit expressions for the Berry phase 7 S in terms of our parametrization of the interaction 
Hamiltonian may easily be obtained. We recall that the unperturbed Hamiltonian and its 
bound and resonant energy eigenfunctions are independent of time. Hence, the time depen- 
dence of the bound and resonant instantaneous energy eigenstates of the perturbed system is 
entirely contained in the matrix K, which is a function of time through the time dependence 
of the external parameters. Therefore, from eq. (50) 

\VrM& Ci(t))) = ]T \^ m (^))(V R K) ms . (63) 

m 

When this expression is substituted in (59), and use is made of the biorthonormality of 
the set of unperturbed bound and resonant eigenfunctions, we get 



Is = 1 



K-^VijK)] -dR. (64) 



The geometric phase is now written as a path integral of the diagonal elements of the product 
of the inverse and the gradient (in parameter space) of the matrix K, which diagonalizes the 
matrix H. This expression for 7 S is equivalent to (59). In the case of a finite number of 
interfering resonant states which may become degenerate, all other bound or resonant energy 
eigenstates being non-degenerate, we may safely assume that H is diagonal except for a square 
diagonal block which mixes the interfering resonances. In this case K is a finite matrix and 
the evaluation of the integrand in the right hand side of (64) involves only the product of two 
finite matrices and no integration over the particle coordinates of the microscopic internal 
components of the system is involved, as would seem to be the case in the evaluation of 
(tp s \V R <p s ). 

Furthermore, from (48) we see that the columns in V^K are the gradients of the instan- 
taneous right eigenvectors of H. Therefore, j s may also be written as 



ls = i J (</> 8 \Vr<I>8) ■ dR (65) 
From (64) we may derive a nice sum rule for the geometric phases of the interfering 



resonante states, 



]T 7s = i / tr [K- 1 (VflK)l • dR. (66) 
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This is a topological invariant, namely, the first Chern class[22], as will be shown below in 
the case of two interfering states. 

Let us consider now the particular case of a system with two resonant states strongly 
mixed by the external interaction, which may become degenerate by a small variation of the 
external parameters while all other bound or resonant states remain non-degenerate. In this 
case the system is in a domain in parameter space which contains one and only one closed line 
of singularites of K which is topologically equivalent to the diabolical circle corresponding to 
an accidental degeneracy of the two interfering resonant states. Then, we may safely assume 
that H is diagonal except for a 2 x 2 block H2 X 2- The matrix K which diagonalizes H2 X 2 is 
given in (36) 

The matrices K -1 and VrK are readily obtained from (14), and (36) 




(67) 



and 



VK = 7?S {-7Te K + ^ Ml ) ( A - *3 f ) + 2T M2? 



(68) 




where the matrices Mi, Mi and M3 are 




(69) 




(70) 



and 




(71) 



Then, a straightforward calculation gives 




(72) 



and 
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1 

72 = -o 



2 7c Te (e - ??) 





(73) 



These expressions are very similar to the well known results obtained for the geometric 
phase of bound states [1]. An obvious difference is that the geometric phase of resonant states 
is complex since e and r\ are complex functions of the parameters R and T. There is another 
important but less apparent difference: In the case of an accidental degeneracy of resonances 
(r / 0), the denominator in the right hand side of (72) and (73) vanishes on the continuous 
line of singularities we have called the diabolical circle, and not at one isolated point as is 
the case for bound states (r = 0). It follows that two kinds of non-trivial, topologically 
inequivalent closed paths are possible, as shown in Fig. 3. First, those paths which surround 
or encircle the diabolical circle but are not linked to it. Second, the closed paths which are 
linked to the diabolical circle. Paths of the first kind are clearly analogous to the non-trivial 
paths that go around the diabolical point while paths of the second kind have no analogue 
in accidental degeneracies of bound states. 

For paths of the first kind it is always possible to find a surface £ which spans the closed 
path c and does not cross the diabolical circle. Then, using Stokes theorem we may write 
the geometric phase 7 S as a surface integral. Computing directly from (72) and (73), we get 



upwards, then Si is above c and £2 is below c. 

This is, of course, the same result as would have been obtained from the general expression 
(62) and our parametrizacion of the perturbation term in the Hamiltonian, no summation 
over intermediate states occurs in (74) since in the simple case of only two interfering resonant 
states the summation in (62) has only one term. 

Adding 71 and 72, the sum rule is written as a surface integral 




(74) 



where s = 1,2. Since 72 changes into 71 when Si and £2 ar e exchanged and the sign of dS 
is changed, the normals for £2 and £1 should be oppositely oriented. If we say that F points 




(75) 
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where £ is a closed surface with the diabolical circle completely contained in its interior, dS 
is the surface element normal to S. 

The integral in (75) is easily computed when £ is a sphere with radius R > ^T, the result 

is 



7i +72 



-2vr 



(76) 



It is now easy to show that the resonance degeneracy produces a continuous distribution 
of singularities on the diabolical circle. To this end, we convert the surface integral (75) to a 
volume integral using Gauss theorem. Then, 



7i + 72 



v 



(i?-^f 



dV, 



(77) 



where V is the volume inside J2 an d bounded by it. The term in round brackets under the 
integration sign vanishes when e / 0. Therefore, the non-vanishing value of 71 + 72 implies 
the occurrence of 5- function singularities of the integrand on those points where e vanishes. 
Hence, 



R-i\T 



(R-i\V 



R 2 



5 (cos 9) 



(78) 



the factor R 2 multiplying the delta function is needed to reproduce the value 2ir of 71 + 72 • 

We may say that, instead of having the fictitious magnetic charge on the monopole singu- 
larity characteristic of the accidental degeneracy of bound states, in the case of an accidental 
degeneracy of resonant states the fictitious magnetic charge is evenly and continuously dis- 
tributed on the closed line of singularities we have called the diabolical circle. 

Now, let us consider the difference between geometric phases of resonant and bound states. 
With this purpose in mind we rewrite (72) and (73) as 



1 r(fxRj-dR 1 rri(fxRj-dR 
71,2 = "2 J c r (e 2 - r? 2 ) ± 2 J c re(e 2 -?? 2 ) 



(79) 



These expressions take a simple and transparent form when we change from cartesian 
coordinates (X,Y,Z), with OZ parallel to T, to spherical coordinates (R,0,ip) in parameter 
space. Then, (79) becomes 
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1 ,- I , i i?cos 8 — ihT) dp 

71,2 = -- dp T ^ / 1 o > — (80) 

* Jc * Jc . / R2 _ 1 r2 _ ,T P r*r>a 



i? 2 -±r 2 -;ri?cos# 

the path c is specified when i? and are given as functions of p. 
The sum rule for the geometric phases takes the form 



7i +72 



■J dtp (81) 



This result is valid for all paths c, it is also valid for bound or resonant states. 
For closed paths of the first kind which go once around the diabolical circle 

7i + 72 = — 27r c* 7 ) of first kind (82) 

This result is valid for all paths , independently of the shape of . It is of course, 
the same result we obtained using the surface integral representation. 

Still in the case of closed paths of the first kind, we may rearrange (80) as 

71,2 = ~\ I [1 T cos 9] dp ± A 7 (83) 

where 



i r [i\T -2Rcos9) + [R- JR 2 -\T 2 -iTRcosO) 
A 7 = i-Y / } 2 I i ' , dp (84) 

4 JcO ^ R 2 _ l F 2 _ iRT CQS g R _ ^ R 2 _ l F 2 _ iRT cos d 

The first term in the right hand side of (83) is the well known expression for the Berry 
phase of two interfering bound states adiabatically transported in parameter space around 
an accidental degeneracy. The second term gives the difference of the actual Berry phase of 
the resonant state and the phase of a bound state, 

7i,2 S = 7i" ± A 7 (c^), of the first kind (85) 

as shown in (84), A7 is proportional to V and vanishes for vanishing T. 

Therefore, the Berry phase acquired by a resonant state when it is transported in pa- 
rameter space around the diabolical circle in a path not linked to it is equal to the sum of 
the real geometric phase a bound state would have acquired if transported around the same 
path, plus a complex correction term characteristic of resonant states. 
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In the case of closed paths of the second kind, that is, those paths which are linked to the 
diabolical circle, there is no surface £ which spans the closed path c^ 11 ^ without crossing the 
diabolical circle, see Fig 3. Therefore we may not use Stokes theorem to convert the path 
integral into a surface integral. However, we may still compute the geometric phase from the 
path integral, 

For closed paths which are linked to the diabolical circle, the sum rule gives zero 

7i + 72 = — / dip = 0, c*- 77 ) of second kind (86) 

since, in this case, the angle tp starts out at some initial value ipo and, as the system traces 
the path c^ n \ it oscillates between a minimum and maximum values and finally ends at the 
same initial value po. 

The difference of 71 and 72 may be obtained from (80), 

71,2 = ±A 7 ( C ( /7 )) (87) 

where A7(c (7/ )) is given by an expression similar to ( 84). There is no analogue to this case 
in bound states. 



5 Validity of the Adiabatic Approximation 

The present calculation of the Berry phase of a resonant energy eigenstate is based upon the 
adiabatic approximation. In general, there seems to be some incompatibility between the 
decay of the system, that is, the vanishing of the signal and adiabaticity, i.e. slow motion. A 
rough estimation of the validity of the adiabatic approximation in this case may be obtained 
from the criterion given in A. Messiah[20], 



maximum angular velocity of \(p s 



minimum Bohr frecuency of \ip s ) 
which, in our notation means 



< 1 



(88) 



ftmax{|E m ^|^f)|} 



min{|£ s - £ m \} 



Now, recalling the relation 



(89) 



Ps 



, d<Prr. 
' dt 



dH x 
~dT 



<Pr, 



(90) 
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and calling 



AE = min j|f s - £ m \] 



(91) 



the condition (89) becomes 



/imax ||^s 



dHi 
dt 



<Pr, 



(AE) 

If T is some typical time of the driving Hamiltonian, such that 



(92) 









1 


($8 


dt 




~ T 



«|ffl|^m)| (93) 

T should not be so long that the signal cannot be meassured. Hence, it seems reasonable 
to estimate T from the largest half width of the interfering resonances 



(94) 



Hence, a rough criterion for the validity of the adiabatic approximation for unstable states 
would be 



max{\(0 s \Hi\ (p m }\} 
(AE) 2 



«1, 



(95) 



where max{\((p s \H± \ <p m )\} indicates the maximum value along the path traced by the 
system in parameter space. 



6 Results and conclusions 

The purpose of the foregoing has been to discuss the adiabatic evolution of an open quantum 
system in a state which is a superposition of resonant states and evolves irreversibly due to 
the spontaneous decay of the unstable states. More specifically, we studied the geometric 
phase acquired by the resonant states when they are adiabatically transported in parameter 
space by the mixing interaction around a degeneracy of resonances . 

In the case of two resonant states mixed by a Hermitian interaction we find two kinds of 
accidental degeneracies which may be characterized by the number and lenght of the cycles of 
instantaneous energy eigenfunctions at the degeneracy. In the first case there are two linearly 
independent eigenfunctions belonging to the same repeated energy eigenvalue, that is, two 
cycles of lenght one. In the second case there is only one degenerate resonant eigenstate 
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and one generalized resonant eigenstate belonging to the same degenerate (repeated) energy 
eigenvalue, i.e. one cycle of lenght two. 

Accidental degeneracies of the first kind, or first rank, give rise to one ' monopole' point 
singularity at a diabolical point in parameter space, as in the case of degeneracies of negative 
energy bound states. In the present case the degenerate states are bound states of positive 
energy embedded in the continuum [21]. In degeneracies of the second kind, or second rank, 
the fictitious magnetic charge is evenly and continuously distributed on a closed line of sin- 
gularities in parameter space, which is topologically equivalent to a diabolical circle. Only 
second rank degeneracies of two resonances produce a true degenerate resonant state. 

Close to a degeneracy of first rank, the hypersurfaces which represent the real and imagi- 
nary parts of the resonance energies in parameter space are two double cones lying in orthog- 
onal subspaces, with their vertices located at the same point, which for this reason might be 
called a double diabolical point. 

When the degeneracy is of the second rank, the topology of the energy surfaces is different 
from that at a crossing of bound states. The energy surfaces of the two resonant states that 
become degenerate are connected at all points in a circle. Close to the crossing, the energy 
hypersurface has two pieces lying in orthogonal subspaces in parameter space. The surface 
representing the real part of the energy has the shape of a hyperbolic cone of circular cross 
section, or an open sandglass, with its waist at the diabolical circle. The surface of the 
imaginary part of the energy is a sphere with the equator at the diabolical circle. The two 
surfaces touch each other at all points on the diabolical circle. 

In the case of two interfering resonant states, the geometric phase acquired by the resonant 
states when transported around the diabolical circle in a closed path which is not linked to 
it, may be written as the sum of two terms 



would have acquired when transported around a diabolical point in a closed path in the same 
parameter space. The second term is complex, it gives rise to a change of the phase and a 
dilation of the resonant state eigenfunction. Its imaginary part may be positive or negative, 
in consequence, it may produce an amplification or a damping of the wave function which 
may compensate or reinforce the attenuation due to the imaginary part of the dynamical 
phase factor. For long lived, narrow resonances, we may expect A7 to be small compared 
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with rybound^ s i nce j s proportional to the ratio ^ which is roughly proportional to the ratio 
of the width to the real part of the resonant energies. 

When the resonant states are transported in a closed path C 11 which does not go around 
the diabolical circle but is linked to it, the geometric phase they acquire is 

7l ^ (c")=±A 7 (c") (97) 

Since it is not possible to find a continuous surface which spans the closed path C n 
without crossing the diabolical circle, we can not make use of the theorem of Stokes to 
convert the path integral into a surface integral. However, it may readily be computed as a 
path integral from the expression 



f { Z - i\Y) IT x R'\ ■ dR' 
A 7 C" = / \ 2 A > (98) 

v ; JcI1 r y(;?-*if)V+^) 

which is obtained from (80) and (86). As in the previous case, A 7 (c H ^j is complex and 
produces changes of phase and dilations in the resonant state wave function. This case has 
no analogue in bound states. 

The sum of the geometric phases acquired by two interfering resonant states which are 
transported around a degeneracy in a closed path in parameter space is a topological invariant, 
namely the first Chern class [22]. For closed paths of the first kind its value is the "magnetic 
charge" on the diabolical circle, and it vanishes for paths of the second kind. 

In conclusion, we have shown that the Berry phase of resonant states differs in various 
ways from that of bound states. It has some interesting mathematical properties not present 
in the Berry phase of bound states. From the physical point of view it is also interesting 
since it has a new term which produces dilations of the wave function and may give rise to 
observable effects not present in the geometric phase factors of bound states. 
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FIGURE CAPTIONS. 



Fig. 1. Close to a crossing of resonances, the surface which represents the real part of 
the energy difference in parameter space has the shape of a hyperbolic cone of circular cross 
section or diabolo. The two resonances are degenerate at the narrowest cross section or waist 
of the diabolo also called the diabolical circle. 

Fig. 2. The surface which represents the imaginary part of the difference of the two ener- 
gies in parameter space close to a degeneracy is a sphere. The upper and lower hemispheres 
represent the imaginary parts of the two neighbouring energies Ei, E2. The equator coincides 
with the diabolical circle. 

Fig. 3. In the evaluation of the Berry phase of two interfering resonances there are two 
kinds of non-trivial, topologically inequivalent closed paths in parameter space. First, those, 
like which go around the diabolical circle but are not linked to it. Second, those, like 
C^ H \ which turn around the diabolical circle and are linked to it. 
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